Thermal History of Solid "^He Under Oscillation 
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We have studied the thermal history of the resonant frequency of a torsional oscillator containing 
solid ^He. We find that the magnitude of the frequency shift that occurs below ~100 mK is multi- 
valued in the low temperature limit, with the exact value depending on how the state is prepared. 
This result can be qualitatively explained in terms of the motion and pinning of quantized vortices 
within the sample. Several aspects of the data are also consistent with the response of dislocation 
lines to oscillating stress fields imposed on the solid. 

PACS numbers: 67.80.-s, 61.72.Lk, 74.25.Qt 



I. INTRODUCTION 

Over the last several years the torsional os- 
cillator (TO) technique has become a popular 
methodi'^'^'i'^'^'i'^'aaOiii^ to study solid ^He. The 
measurements, which involve monitoring the resonant 
frequency / of a high quality factor torsion pendulum 
filled with ''He, have revealed a "transition" below an 
onset temperature Tq of ^200 mK. The results of 
several control experiments^ indicate that the sudden 
increase in the frequency is due to a superfluid-like 
decoupling of a fraction of the ''He mass. The apparent 
nonclassical rotational inertia fraction (NCRIF), propor- 
tional to the frequency shift 6f, is independent^ of the 
maximum oscillation speed of the TO provided that 
it is smaller than a critical value. Above this critical 
velocity vc the magnitude of NCRIF is attenuated. The 
value of Vc corresponds to several quanta of circulation 
K, suggesting that the important excitations in the 
system are vortices. Anderson has proposed^ a vortex 
liquid model that qualitatively captures a number of 
experimental results. 

In multiple solid ^He samples slowly grown within an 
annulus of radius, r = 5 mm, and width'^'*, t = 0.95 mm, 
a value of vc = 10 ^ms"'^ was found^ (corresponding to 
K = 2'KrvQ Sh/m, where h is Planck's constant and 
m is the bare "^He mass). NCRIF decreased almost lin- 
early with ln[uo] as the speed was increased further. This 
attenuation, according to Ref. [l^, is attributable to the 
nonlinear susceptibility of a vortex liquid phase, which 
consists of an entanglement of many thermally activated 
vorticesi^. The ability of the vortices to move counter to 
the time-dependent superflow (relative to the cell's os- 
cillation) results in the screening of supercurrents. As 
the temperature is lowered the motion and number of 
vortices is reduced so that NCRIF becomes finite. The 
accompanying dissipation^ that peaks in the middle of 
the transition is due to a matching of / and the optimal 
rate at which vortices respond to changes in the veloc- 
ity field. Residual '^He atoms may^^ cling to vortices, 
thereby slowing vortex motion and ultimately enhancing 
the onset temperaturei^iiSiii. 

One prediction of Anderson's modeli^ is that Tq , and 
perhaps NCRIF, for a particular sample will decrease 



as the measurement frequency is lowered. Aoki et al. 
performed^ measurements in a double oscillator that op- 
erated in either an antisymmetric (high /) or a symmetric 
(low /) mode and found that for the same *He sample. 
To « 240 mK at 1173 Hz and Tq « 160 mK at 496 Hz. 
However, NCRIF became independent of frequency be- 
low '^35 mK. In addition, irreversible changes in NCRIF 
below '^45 mK occurred upon variation of the oscillation 
speed between 10 /ims^^ and 800 /ims^^. 

Nonsuperfluid explanations of the TO experiments 
have also been proposed, in which some sort of freezing 
processi^ or glass transitioi*ii takes place. Although the 
former was only discussed qualitatively, the pinning of 
dislocation lines within solid "^He is one possible realiza- 
tion of such a "freezing" transition. The mechanism for 
dislocation pinning is likely related to interactions with 
isotopic impurities. It is known that dislocations'^* and 
''He atom s '^'^° are mobile in solid ''He below ~1 K, but 
the strong coupling between them in the low temperature 
limit can immobilize them both2ii2^. In a TO studyii of 
many solid *He samples we found that the dependence 
of To on the '^He concentration 2:3 is quantitatively con- 
sistent with this impurity-pinning model. 

Recently Day and Beamish observed^^ a dramatic in- 
crease (between 5% and 20%) in the shear modulus C44 
of solid "^He (2:3 ~ 0.3 ppm) below ~250 mK, with the 
temperature dependence of C44 strongly resembling that 
of / in TO experiments. The authors also reported hys- 
teresis in C44 upon adjusting the stress amplitude of the 
measurement at T w 20 mK, similar to what is seen in 
Ref. J,. The increase in the shear modulus is very likely 
due to the stiffening of the dislocation network via the 
immobilization of individual lines. It was suggested by 
the authors of Ref. [2^ that an increase in C44 will result in 
a stronger coupling between the ^He sample and its con- 
tainer and therefore mimic an enhanced mass loading of 
a TO, causing / to decrease. However, any enhancement 
to the overall rigidity of the system actually leads to a 
higher resonant frequency. We have carried out a finite 
element method (FEM) calculation^! employing realistic 
parameters of the TO and solid helium in order to con- 
firm this result. The positive correlation between C44 and 
/, their resemblances in the temperature dependence, and 
the hysteresis common to both, all suggest that either the 
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FIG. 2: Temperature dependence of Sfne and SQhb^^ for the 
same samples presented in Fig. [T] obtained by subtracting the 
temperature dependence of the empty cell. The symbols for 
BC-12/10/06, pre- and post-annealing, are identical to those 
in Fig. [2 



FIG. 1: Temperature dependence of Sf for a BC sample be- 
fore and after annealing, and one CT sample. For all samples 
2:3 ~ 1 ppb. Data was obtained with the beryllium-copper 
TO from Ref . [Tol. Freezing (Tf) and annealing (Ta) temper- 
atures are given for each sample. All data was taken during 
cooling scans in which the temperature was successively low- 
ered in controlled steps. For any one sample the measured 
frequency at T > To is identical for all values of vq- For easy 
comparison we have set 5/ = at T = 175 mK by subtract- 
ing 1071.007 Hz, 1071.011 Hz, 1071.035 Hz, and 1071.898 Hz 
from the data of BC-12/10/06, BC-12/10/06 post-anneahng, 
CT-010507, and the empty TO. Thus, the effective mass load- 
ing for each of these samples is 891 mHz, 887 mHz, and 863 
mHz. 



frequency shift is a direct consequence of the enhanced 
modulus and can be explained without invoking super- 
solidity or that the microscopic mechanism responsible 
for NCRl also affects the elastic properties of the solid 
(e.g., a recent suggestion is given in Ref.|25|). The results 
presented below appear to favor the latter. 

Sec. we present our experimental observations in a 
general manner by simply describing the dependencies 
of the observed frequency shifts and dissipation signals 
on temperature (see Figs. [1] and [2]), oscillation speed-'' 
(see Fig. [3]), and growth method. Two different inter- 
pretations of the data are given in Sees. Illll and IIVI In 
the former we attempt to describe the experiments in 
the context of a supersolid that is permeated by vortices. 
Since vortex lines interact with velocity fields one of the 
important parameters in this model is vq. In Sec. IIVI 
we discuss the consequences of attributing the entire fre- 
quency shift to an increase in the shear modulus of solid 
*He. It is assumed in this picture that the rigidity of 
the solid is strongly affected by the motion of dislocation 
lines under the applied stress a. In anticipation of these 
two interpretations, we have simultaneously plotted the 
data in Fig. [3] as a function of vq and cr. In addition, in 



both Figs, m and Owe have plotted the frequency shift in 
terms of NCRIF and Sc^i/ C44. 



II. EXPERIMENTAL RESULTS 

A recent study^^ of samples grown at constant tem- 
perature and pressure (CT/CP) at a single point on 
the solid-liquid coexistence curve found the onset of the 
frequency shift to be sharp and reproducible, which is 
consistent with the expected formation of single crys- 
tals. In contrast, wide variations in Sf and Tq were 
observed when employing the blocked capillary (BC) 
method, which typically yields polycrystalline samples or 
highly strained crystals due to the pressure drop during 
freezing. In this paper we examine in detail the different 
thermal histories of Sf in BC and CT/CP samples that 
follow from changes in oscillation speed above and below 
To- We have measured the temperature dependence of 
Sf in 30 samples, with the bulk of our work being carried 
out on "isotopically pure" (0:3 « 1 ppb of '^Hc) samples at 
1 /.tms^^ < vo < 100 /ims~^. To overlap with Ref.Q we 
also studied commercially pure (xa ~ 0.3 ppm) samples 
up to speeds of 880 ^ms~^. 

Figure [1] shows the temperature dependence of Sf for 
a BC sample before and after annealing, compared with 
that for a sample grown at CT. We find that the temper- 
ature dependence of Sf for each sample is reproducible 
in warming and cooling scans at low oscillation speeds, 
i.e., on the order of one micron per second. Although the 
differences between BC and CT/CP samples are appar- 
ent at low speeds, the temperature dependence obtained 
at Vq > 50 iiuis~^ is very similar in all samples studied 
(including 0.3 ppm samples). 

After subtracting the temperature dependence of the 
empty cell we obtain the frequency shift and dissipation 
due to solid ^He, which we respectively denote by SfHe 
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FIG. 3: Dependence of the low temperature value of Sfne 
on oscillation speed and applied stress. The inertial stress 
is approximated by o" « Trpfvo/2, for a given helium density 
p. Values of fne were extracted from cooling scans that be- 
gan well above To- Samples are labeled according to growth 
method and date. The data are plotted versus a linear ab- 
scissa in the inset. 



and SQhc^^- These results are plotted in Fig. O There 
is considerable dissipation in rapidly grown BC samples 
even at the lowest temperatures. After annealing sam- 
ple BC-12/10/06 the high temperature tail of Sfne (and 
thus To), the low temperature limiting value of (5/ffe, and 
the magnitude of SQHe~^ are all reduced. Thus it ap- 
pears that annealing affects the temperature dependence 
of Sfne much in the same way as driving the TO at high 
speed (see Fig. [5]). One difference, however, is that the 
results from annealing are permanent. 

In Fig. [3] we have plotted the low temperature value 
of Sfne measured at different vq, which again demon- 
strates that the data converge at high speeds. One in- 
teresting difference from previous TO studie9^i^>^iii2^ is 
the lack of saturation in Sfue even well below 10 /ims^^. 
This is most evident in two BC samples (01/25/07 and 
12/10/06) grown from the normal fluid phase, where Sfne 
actually doubles when vq is decreased from 10 ^ms^^ to 
1 ^ms~^. In samples assumed to be of higher quality, 
we find vc ~ 3 ^ms~^ (3.5 /zms^^ corresponds to k = 
h/m). 

The curves in Fig. [3] were compiled from cooling scans 
that began well above the onset temperature. Using this 
procedure we found that Sfne is diminished at higher 
vq for all T < Tq in a completely reproducible man- 
ner, whereas measurements taken during warming are 
history-dependent when obtained at speeds greater than 
a few microns per second. As we will show below, varying 



the oscillation speed below the onset temperature can re- 
sult in irreversible changes in S/hc- Careful examination 
of how the oscillation speed affects Sfne reveals metasta- 
bility at the lowest temperatures, as depicted in Figs. 3] 
and [5] (see Sees. |ni]and|IV|. 

Using the following protocol, the thermal history of 
SfHe was investigated in crystals grown under different 
conditions. First, the temperature dependence of dfne 
was measured while cooling at low speed (vq < 3 /xm s~^). 
At T « 20 mK vq was slowly increased. Final veloci- 
ties were usually ^20 /ims~^, but in some cases much 
greater (e.g., several hundred microns per second). Mul- 
tiple thermal cycles were then performed in succession, 
indicated by the arrows in Fig. U) Complete equilibration 
was purposely avoided while cycling between 30 mK and 
60 mK in order to observe multiple metastable states. 
The measurement culminated with a high velocity cool- 
ing trace that began well above Tq- 

Figures H] and [5] show that for one oscillation speed 
a number of different values of Sfne (« NCRIF and 
6044/044) can be "frozen in" below a characteristic tem- 
perature, T* w 30 mK. This is true regardless of the 
sample growth procedure. However, there are significant 
differences in the time evolution (following temperature 
steps) and the velocity dependence of Sfne between sam- 
ples formed by the BC method and crystals grown at 
CT/CP. Similar behavior is seen in 0.3 ppm samples, but 
with r* « 45 mK. For all samples, there is no observ- 
able decay from any of the metastable states of Sfne on 
a timescale of days as long as the temperature is main- 
tained below T*. Further discussion of the hysteresis is 
deferred to Sees. Illll and IIVI below. 



III. THE VORTEX LIQUID MODEL 

As stated in Sec. U the TO experiments to date are 
qualitatively consistent with many features of Anderson's 
vortex liquid model. The present work, however, demon- 
strates that below 60 mK the system does not act like 
a free vortex liquid as originally proposed. In fact, be- 
low 30 mK it exhibits severe pinning. With this picture 
in mind we can explain the temperature dependence of 
NCRIF (oc 6fHe) in Fig. [21 the velocity dependence in 
Fig. [21 and the history dependence in Figs. [Hand O 

First we consider the temperature dependence. As the 
system is cooled through the transition some number of 
free vortices are produced, which is proportional to vq. 
When T « To the mobility of the existing vortices di- 
minishes and it is unfavorable to create additional vortex 
lines. The inability of vortices to screen supercurrents 
leads to a sizeable NCRIF in the low temperature limit. 
When the system is driven at high speeds, not only is 
NCRIF smaller but the high temperature tail of NCRl is 
reduced. This suggests that vortex pinning is very weak 
at T > 60 mK. It is also likely that vortices are more 
easily excited as To is approached. 

The velocity dependence in Fig. [3] reflects the distri- 
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FIG. 4: Metastability of Sfne (plotted as NCRIF and 
JC44/C44) for a BC sample. The values of 5C44/C44 were cal- 
culated based on the assumption that the solid sample was 
isotropic, i.e., highly polycrystalline^^ due to BC method that 
was employed. Even larger values would be obtained if the 
sample was assumed to be a single crystal. Decay of the sig- 
nal occurs for T > 30 mK. Cycling above and below T* — 
30 mK results in successively smaller values until equilibrium 
(i.e., that obtained when cooling from above To) is reached. 
The speeds of vo = 1.3 /ims~^ and 21 /^ms"^ correspond to 
(7 = 2.0 mPa and 32 mPa. 



bution in the strength of individual vortex pinning sites 
at low temperature. Annealing might effectively remove 
the weakest pinning centers so that the low temperature 
limiting value of NCRIF is only decreased at small wq. 
At greater speeds these weak regions are irrelevant since 
there is a plethora of free vortices, hence the convergence 
in Figs. [2] and [3] at > 50 /ims^^. 

In the context of superconductivity or magnetism (or 
superfluidity), the system is prepared differently if (veloc- 
ity) field-cooling or zero field-cooling are employed. The 
hysteresis shown in Figs. [3] and [5] is a consequence of what 
the authors of Ref. have called vortex glass behavior. 
If the state of the system is prepared at low vq there is 
only a "small" number of vortices frozen into the solid. 
When vq is gradually increased while T < T*, NCRIF 
becomes either unstable (see Fig. 01) or metastable (see 
Fig. [5]). In the former case, NCRIF decays until it reaches 
a metastable state. For T > T*, NCRIF is reduced by 
the enhanced mobility of vortices. 

The low temperature metastability seems to depend 
on the difference between the initial NCRIF prepared 
at low speed and the "equilibrium" value obtained upon 
cooling at the higher speed. For instance, the difference is 
NCRIF = 0.9% for BC-12/10/06 (see Fig.© and NCRIF 
= 0.08% for CP-11/21/06 (see Fig.^, resulting in a much 
more stable value for the latter. The measured value 
for BC-12/10/06 immediately drops by about 40% of its 
magnitude following the increase in vq, while there is no 



FIG. 5: Metastability of NCRIF and SC44/C44 for a CP sam- 
ple. The sample was assumed to be a single crystal with its 
c-axis oriented horizontally. If we instead assume that the 
c-axis is oriented vertically (or that the sample is isotropic) 
the low temperature limiting value of 5C44/C44 would be 60% 
(30%). Labels A through I denote changes in the decay rate at 
several different temperatures. The speeds of vo ~ 2.2 /ims~^ 
and 22 /ims~^ correspond to cr = 3.4 mPa and 34 mPa. 



change for CP-11/21/06. The behavior of CP-11/21/06 
is similar to that in Ref. 0, which was compared with the 
Meissner effect (i.e., a robust NCRIF due to the inability 
of vortices to enter into the solid). We note that the 
metastability up to 800 /j,ms~^ in Ref.Qmay be related to 
the small NCRIF « 0.1% of the sample. In comparison, 
upon raising vq to 880 /ims~^ we observe, for a 0.3 ppm 
sample with NCRIF > 1% {Sfne > 8 mHz), an abrupt 
drop in the signal that is followed by a gradual decay. 
In fact, in this extreme limit of vq, the decay rate is 
essentially insensitive to temperature steps (warming and 
cooling) between 20 mK and 45 mK. Only upon warming 
above 60 mK does the decay rate increase further. 

In 1 ppb samples, the stability is reduced as the tem- 
perature is raised above 30 mK. Fragments of data from 
scans of CP-11/21/06 are displayed in Fig. [6] The data 
at any T < 60 mK cannot be fit satisfactorily using a 
simple mathematical formula. For example, the decay 
at 40 mK in Fig. [Ifa) (between points A and B) begins 
immediately, but slows down abruptly after only 1 h. 
In Fig. E^b), NCRIF barely changes in the first hour at 
the same temperature (data prior to point G). Irregular 
rates are even more apparent at 50 niK. The decay slows 
down dramatically when NCRIF crosses (e.g., points E 
and H) from the region above the low velocity cooling 
trace to the region below it. This is most obvious in 
the Fig. [Bfb), where the sample was warmed to 50 mK 
rapidly to result in a larger initial NCRIF. The very dif- 
ferent decay rates in these two regions, and the identical 
temperature dependence in warming and cooling scans, 
both suggest that there is physical significance to the re- 
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FIG. 6; Segments of raw data used to compile the "cycling" 
and "warming" curves in Fig. [5] Decay rates above (D to 
E and G to H) the low velocity trace are greater than those 
below (E to F and H to I). The ordinate axis is given in terms 
of NCRIF for easy comparison with Fig. [S] 



producible data obtained in the low velocity limit. Simi- 
lar behavior could not be observed in BC samples because 
of higher decay rates, i.e., it was impossible to exceed the 
low velocity trace upon warming (see Fig. 3]). 

Interestingly, the decay in CP-11/21/06 at 60 mK 
(data following point I) is different from that at 40 mK 
and 50 mK in that it can be well fit to an exponen- 
tial form, which yields a time constant of 2.25 h. The 
smooth decay across the low velocity trace may indicate 
that the small tail of NCRIF for T > 60 mK does not de- 
note the same "boundary" as seen at lower temperatures. 
This observation is complemented by the enhanced de- 
cay rate at 60 mK that we find in 0.3 ppm samples as 
well as the ss-independent specific heat peak that has 
been observedlS in the same vicinity, and suggests that 
there is some inherent transition between the "free" and 
"pinned" regimes. 

Another intriguing result is the nearly equal spacing 
between each metastable NCRIF among the four sam- 
ples that we studied in detail (e.g., see Figs. [4] and [5]). It 
is clear that the system prefers to possess specific values 
of NCRIF in the low temperature limit. For example, 
counter to the slow decay upon warming a solid sample, 
as the temperature is lowered we find NCRIF to "jump" 
up to a steady value (see the data in Fig.[5]following point 
F) . If we naively estimate the frequency shift due to the 
presence of a single vortex line using the expression^. 



TABLE I: Approximate stress amplitudes corresponding to 
vc in various cells. Values are estimated by ignoring the ef- 
fects from the top and bottom surfaces of the cell and equat- 
ing two different expressions for the torque: la — rAa. Here, 
a is the angular acceleration and A is the area at the rim 
of the cell. The stress exerted in an annular geometry dif- 
fers from that in a cylindrical one by the approximate factor, 

{2tf)Ann/{rf)cyl. 



Cell 


./ [Hz] 


vc [fims ^] 


ac [mPa] 


Cylinder" 


1072 


< 3.5 


< 5.4 


Cylinder' 


1173 


15 


12 


Cylinder* 


496 


15 


28 


Annulus'^ 


912 


3-40 


1.6-22 


Anrmlus'' 


912 


10 


5.4 


Annulus*^ 


874 


40 


3.3 



"Present set of measurements. 
''Data from Ref. 
"^Data from Ref. 
''Data from Ref. 
■^Data from Ref. 



= 7rpr4NCRIF/ln[r/ao], we calculate a shift cor- 
responding to 25% of the observed NCRIF for a vortex 
core radius^- of ao =0.1 nm. This value is of the same 
order of magnitude to what we observe (e.g., from Fig. [5] 
we get 0.04/0.33 = 11% for CP-11/21/07). 

Although our experimental results can be qualitatively 
interpreted within the framework of Anderson's vor- 
tex liquid model, very few quantitative comparisons are 
possible^. One of the present experimental results that 
is difficult to understand in the vortex picture is depicted 
in Fig. [3l In two samples the observed frequency shift 
continues to increase when vq is reduced to 1 /ims~^, 
such that even the peak to peak amplitude of 2 ^ms~^ 
is less than that which corresponds to a single quantum 
of circulation (^^3.5 /ims~^). 



IV. ANOMALOUS ELASTICITY OF SOLID ^HE 

Due to the similarities discussed in Sec. [J between the 
observed frequency shifts and the results of Ref. [2^ we 
are compelled to interpret the data from Sec. |TT]in terms 
of the response of the dislocation network to oscillating 
stress fields imposed on the solid. 

Here we report our preliminary estimates of changes in 
the shear modulus of solid ^He. By calculating^'* the reso- 
nant frequency of the torsion mode for different C44 values 
we extracted the percentage change in the modulus that 
corresponds to a small frequency shift of the TO. The 
change is given by (5c44/c44 = {dlnf / d\iiC44y^{SfHe/ f), 
where the derivative term comes from the FEM calcu- 
lation and frequency shift comes directly from our TO 
measurements. The results of our calculations are shown 
in Figs, m and [HI Since the low temperature change in the 
resonant frequency is very small (parts per million), it 
scales linearly with the change in the shear modulus. As 
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a result the inferred temperature dependence of 5C44/C44 
is identical to Sfne ■ It is therefore natural that the qual- 
itative features of the data in Ref. [13 be very similar to 
that from TO studies. 

However, the magnitude of the relative change in C44 
appears to be different. In order to account for the ob- 
served 6fHe in TO measurements the required 6044/0^4 
values among the samples studied are between two and 
five times larger than those seen in Ref. Moreover, 
the changes are greater than the theoretical upper limit 
of 30% that can result from the pinning of dislocation 
lincs'^^. We have carried out similar calculations for some 
of the TO's in the literature^^i '^i ' Ti . ?- ! and found that even 
larger values of (5c44/c44 are necessary to account for the 
observed 6fHe- For instance, typical NCRIF values in 
Refs. i2u4432i all translate into an increase in C44 by ap- 
proximately a factor of 10. For the double oscillator from 
Refl3 we find C44 to increase by factors of 5 and 10 in 
the anti-symmetric and symmetric modes, respectively. 
This implies that the effect is very dependent on fre- 
quency, contrary to the findings of Day and Beamish. 
Furthermore, it is difficult to apply the dislocation pin- 
ning picture to the experiments in porous Vycor glass^, 
in which the dimensions of helium (i.e., the pore size) 
are smaller than typical dislocation loop lengths in bulk 
crystalsi^i2ii2^. This, and the fact that the crystalline 
phased of the confined ^He is body-centered-cubic rather 
than hexagonal-close-packed, would lead one to expect 
very different behavior in the response to torsional oscil- 
lations, counter to what is observed. 

Day and Beamish found the measured shear modulus 
to be a function of stress amplitude. Below a critical 
stress ac the magnitude of C44 saturates at a maximum 
value, as does S/hc in some samples (see Fig. [3]). This 
critical value, ac = 300 mPa, is smaller than the esti- 
mated breakaway stress of 4 Pa for a ^He-^He separation 
of 5 /im (a typical loop Iengtbi^i2ii2^) along a dislocation, 
but closely matches the stress at which S/hc tends to zero 
in Fig. [3] Such strong nonlinearity between the drive 
and response, also seen in other low frequency acoustic 
measurements--, occurs at "critical" speeds that are 100 
to 1000 times less (e.g., ~50 nms^^ in Ref. (231) than in 
TO studies^i^i^iii^. It was concluded in Ref. [H that ac is 
a more fundamental quantity than vc due to it being in- 
dependent of the measurement frequency. However, this 
is inconsistent with Ref. 7 and implies that the two types 



of measurements, although clearly related on some level, 
are distinct. 



V. CONCLUSIONS 

We have studied a number of solid *He samples grown 
within a torsional oscillator at constant temperature and 
pressure, as well as with the blocked capillary method, 
and found that the resonant frequency of the system pos- 
sesses a very strong thermal history. The magnitude 
of the apparent NCRIF measured in the low tempera- 
ture limit is reproducible when obtained upon cooling 
the sample from temperatures well above the transition. 
Modulation of the oscillation speed below the onset tem- 
perature reveals the existence of many metastable states. 

We find that these properties are qualitatively consis- 
tent with two different interpretations: the response of 
vortices to velocity fields and/or the response of dislo- 
cation lines within the solid to oscillating stress fields. 
The mobility of either of these entities becomes very lim- 
ited with decreasing temperature, resulting in hysteretic 
behavior below 60 mK. Extreme pinning takes place be- 
low 30 mK (1 ppb) or 45 mK (0.3 ppm), such that their 
motion is essentially frozen out. 

At present there remain discrepancies between the 
experiments and either of the two interpretations. In 
the vortex liquid picture, the most notable problem is 
the very low critical velocity found in this study. In 
regard to the elastic properties of solid '*He, we find 
that the majority of the frequency shifts observed are 
appreciably larger than upper limit expected from the 
pinning of dislocations. 
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